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Abstract 

The dynamical model developed in [Phys. Rev. C 54, 2660 (1996)] has been 



£> ' applied to investigate the pion electroproduction reactions on the nucleon. It 

in ' 

is found that the model can describe to a very large extent the recent data of 
p(e,e'ir°) reaction from Jefferson Laboratory and MIT-Bates. The extracted 

o" 

magnetic dipole(Ml), electric dipole(E2), and Coulomb(C2) strengths of the 
7JV — * A transition are presented. It is found that the C2/M1 ratio drops 
significantly with Q 2 and reaches about -13% at Q 2 =4 (GeV/c) 2 , while the 
E2/M1 ratio remains close to the value ~ —3% at the Q 2 = photon point. 
^ ' The determined Ml transition form factor drops faster than the usual dipole 



form factor of the proton. We also find that the non-resonant interactions 

can dress the 7./V — > A vertex to enhance strongly its strength at low Q 2 , but 

much less at high Q 2 . Predictions are presented for future experimental tests. 

Possible developments of the model are discussed. 
PACS number(s): 13.75.Gx,21.45.+v,24.10-i,25.20.Lj 
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I. INTRODUCTION 



One of the purposes of investigating the nucleon resonances(iV*) is to understand the 
non-perturbative dynamics of Quantum Chromodynamics(QCD). One possible approach to 
realize this is to compare the predictions of QCD-inspired models with the resonance param- 
eters which can be extracted from the data of pion photoproduction and electroproduction 
reactions. In recent years, precise data including polarization observables have been obtained 
in the A region for pion photoproduction at LEGS [|TJ and Mainz M, and for pion electro- 
production at Thomas Jefferson National Accelerator Facility(JLab) MM, MIT-Bates |5| 
and NIKHEF f| . These data now allow us to investigate more precisely the electromagnetic 
excitation of the A resonance. 

In Ref. we have developed a dynamical model(henceforth called the SL model) to 
extract the magnetic dipole(Ml) and electric quadrapole(E2) strengths of the 'jN — > A 
transition from the pion photoproduction data. The precise polarization data from LEGS 
and Mainz were essential in our analysis. In this paper, we report on the progress we have 
made in extending the SL model to investigate the pion electroproduction reactions in the A 
excitation region. We will make use of the recent data from JLab and MIT-Bates to explore 
the Q 2 -dependence of the jN — > A transition and make predictions for future experimental 
tests. 

The dynamical content of the SL model has been given in detail in Ref. . The essential 
feature of the model is to have a consistent description of both the ttN scattering and the 
electromagnetic production of pions. This is achieved by using a unitary transformation 
method to derive an effective Hamiltonian defined in the subspace ttN © 7iV © A from the 
interaction Lagrangians for N, A, n, p, uj and photon fields. The resulting model has given 
a fairly successful description of the very extensive data for pion photoproduction. The 
extension of the SL model to investigate pion electroproduction is straightforward. The 
formulae needed for calculating the current matrix elements of pion electroproductions are 
identical to that given in Ref. 0] except that a form factor must be included at each photon 
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vertex. Therefore no detailed presentation of our model will be repeated here. Similarly, 
we will not give detailed formulae for calculating the electroproduction cross sections since 
they are well documented |^I2|. 

The SL model is one of the dynamical models developed [|T|,|T3|-f2TH in recent years. Com- 
pared with other approaches based on the tree-diagrams of effective Lagrangians 



or 



dispersion-relations [p5|-p^|, the main objective of a dynamical approach is to separate the 
reaction mechanisms from the excitation of the internal structure of the hadrons involved. 
Within the SL model, this has been achieved by applying the well-established reaction theory 



within the Hamiltonian formulation (see, for example, Ref. ||28|| ). In particular, the off-shell 
non-resonant contributions to the 'yN — > A form factors can be calculated explicitly in a 
dynamical approach. Only when such non-resonant contributions are separated, the deter- 
mined "bare" 'yN — > A form factors can be compared with the predictions from hadron 
models. Within the SL model, this was explored in detail and provided a dynamical in- 
terpretation of the long-standing discrepancy between the empirically determined magnetic 
Ml strength of the — > A transition and the predictions from constituent quark mod- 
els. In this work, we further explore this problem utilizing the Q 2 -dependence accessible to 
electroproduction reactions. Furthermore, the Coulomb(Scalar) component C2(S2) of the 
7iV — > A form factor will be determined. 

In section II, we briefly review the essential ingredients of the SL model and define 
various form factors which are needed to describe pion electroproduction reactions. With 
the Mainz data 0], we have slightly refined our model at the Q 2 = photon point. This 
will be reported in section III. The electroproduction results are presented and compared 
with the data in section IV. In section V, we give a summary and discuss possible future 
developments. 
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II. THE SL MODEL 



Within the SL model, the pion photoproduction and electro-production reactions are 
described in terms of photon and hadron degrees of freedom. The starting Hamiltonian is 

H = H + Hj, (1) 

with 

Hi = ^2 Fmb+->b', (2) 

M,B,B' 

where Hq is the free Hamiltonian and Tmb^b 1 describes the absorption and emission of 
a meson(M) by a baryon(B). In the SL model, such a Hamiltonian is obtained from 
phenomenological Lagrangians for N, A, 7T, p, to and photon fields. In a more microscopic 
approach, this Hamiltonian can be defined in terms of a hadron model, as attempted, for 
example, in Ref. . 

It is a non-trivial many body problem to calculate nN scattering and 'yN — > ttN reac- 
tion amplitudes from the Hamiltonian Eq. ([I]). To obtain a manageable reaction model, 
a unitary transformation method |]7|j30| is used up to second order in Hi to derive an ef- 
fective Hamiltonian. The essential idea of the employed unitary transformation method is 
to eliminate the unphysical vertex interactions MB — > B' with itim + wis < m B' from the 
Hamiltonian and absorb their effects into MB — > M'B' two-body interactions. In the SL 
model, the resulting effective Hamiltonian is defined in a subspace spanned by the 7rN, 'yN 
and A states and has the following form 

Heff — H + v w n + v 11T + T nN ^A + r 7 jv^A, (3) 

where v^n is a non-resonant tcN potential, and t> 77r describes the non-resonant 7iV <-> ttN 
transition. The A excitation is described by the vertex interactions T^^a for the 7iV <-» 
A transition and T^n^a for the ttN <-> A transition. The vertex interaction r 7 A^A is 
illustrated in Fig. 1. The non-resonant v 17T consists of the usual pseudo- vector Born terms, 
p and uj exchanges, and the crossed A term, as illustrated in Fig. 2 (the non-resonant term 



due to an intermediate anti-A state was found to be very weak and can be neglected). Most 
of the dynamical models have the above form of the Hamiltonian. However, the SL model 
has an important feature that the deduced effective Hamiltonian H e ff is energy independent 
and hermitian. Hence, the unitarity of the resulting reaction amplitudes is trivially satisfied. 
Furthermore, the non- resonant interactions v 77r and v n pf are derived from the same Hj of Eq. 
(H), and hence the 11N and reactions can be described consistently. Such a consistency is 
lost if v n N is either constructed purely phenomenologically as done in Refs. [|l"3HT5|l or taken 



from a different theoretical construction as done in Refs. |[16|j21| . This consistency is essential 
in interpreting the extracted 7iV — > A form factors since the non-resonant interactions tViv 
and dress the 'jN — > A vertex. As discussed in Refs. 0,0, only the dressed — > A 



transition can be identified with the data. The importance of the non-resonant effects on 



the jN — > A transition is also stressed recently in Ref. |H . 

From the effective Hamiltonian Eq. (0), it is straightforward to derive a set of coupled 
equations for nN and reactions. The resulting pion photoproduction amplitude can be 
written as 

=< 71N I e- J I 7N >, (4) 

where J is the current operator and e is the photon polarization vector. It can be decomposed 
into two parts 

The non-resonant amplitude t 77r is calculated from t> 77r by 

tyx(E) = f 77r + t 7rN (E)G 7rN (E)v 17T , (6) 
where the irN free propagator is defined by 

° ME) = E-E N {k) l -EM +l e- (7) 
The amplitude t^N i n Eq. (^|) is calculated from the non-resonant irN interaction v^n by 
solving the following equation 



tnN(E) = V-irN + V^G^{E)t^{E) . (8) 

The dressed vertices in Eq. @ are defined by 

T 7 a^a = T 7 jv-»A + r vr Ar^AG' 7r Ar(-E)f 77r , (9) 
Ta^ttTV — [1 + ^7rAr(-E')G vr Ar(i?)]rA^7rAr- (10) 

In Eq.(§J), we also have defined 

[1 + G nN (E)t nN (E)]. 
The A self-energy in Eq. (|5|) is then calculated from 

Ea(-E) = r 7rA r^AG' 7r iv(-E)rA- > 7rAr. (11) 

As seen in the above equations, an important consequence of the dynamical model is that 
the influence of the non-resonant mechanisms on the resonance properties can be identified 
and calculated explicitly. The resonance position of the amplitude defined by Eq. ([|) 
is shifted from the bare mass ira by the self-energy £ a (-£/)• The bare vertex r 7 jv-*A is 
modified by the non- resonant interaction t> 77r to give the dressed vertex T^ N ^^, as defined 
by Eq. (^|). In the SL model, it was found that the extracted Ml strength of the bare vertex 
r 7 A^A is very close to the values predicted by the constituent quark models [p9|-|4T|], while 



the empirical values given by the Particle Data Group(PDG) []32| can only be identified with 
the dressed vertex T^ N ^^- 

The above equations can be solved for arbitrary photon four-momentum Q 2 = —q 2 > 0. 
For investigating the electroproduction reactions, we only need to define a form factor at 
each photon vertex in Figs. 1-2. For the non-resonant interactions(Fig. 2), we follow the 



previous work flQjlTlj . The usual electromagnetic nucleon form factors(given explicitly in 



the Appendix A of Ref. Jl0|) are used in evaluating the direct and crossed nucleon terms. 



To make sure that the non- resonant term t> 77r is gauge invariant, we set 

F A (q 2 ) = F^ 7r (q 2 )=Fl / ( q 2 ), (12) 



where Fa{q 2 ) is the form factor for the contact term, F^ niT (q 2 ) is the pion form factor for the 
pion-exchange term, and FY{q 2 ) is the nucleon isovector form factor (also given explicitly 
in the appendix A of Ref. flOfl). The form factors for the vector meson-exchange terms are 
chosen to be 

gvnjiq 2 ) = gvfnl i 1 ~ q 2 / m v), (13) 

where m v is the vector meson mass and the coupling constants gvvy for V = p,u are deduced 
from the V — > -yn decay widths and are given in Ref. 0. The prescriptions Eqs. (p~2|) - (p~3|) 
have been commonly used in most of the previous investigations such as that in Refs. |Kj[lI|. 
Undoubtly, this is an unsatisfactory aspect of this work. On the other hand, " dynamically" 
sound progress in solving the gauge invariance problem cannot be made unless a microscopic 
theory of hadron structure is implemented consistently into our model. This is beyond the 
scope of this work. 

For the 'jN — > A form factors, we extend Eq. (4.15) of Ref. to write in the rest frame 
of the A 



< A | T 7 at^a | q > 



E N (q)+m N 1 3(m A + mAr) 

-*3 



(2tt) 3 / 2 ^ 2E N (q) ^Am N (E N (q)+m N 
x [iG M (q 2 )Sx q-e + G E (q 2 )(S-ea-q + S-qa-e) 
, G c (q 2 " 



-S-qcj-qe ], (14) 



where e = y An/137, q = (u, q) is the photon four-momentum, and e = (eo, e) is the photon 
polarization vector. The transition operators S and T are defined by the reduced matrix el- 
ement < AjlS'lliV >=< A||T||iV >= 2 in Edmonds' convention ||33|| . The parameterizations 



of the form factors GM{q 2 ), GE^q 2 ) and Gc(q 2 ) will be specified in section IV. 

With the form factors defined in Eqs.(12)-(14), both the non-resonant term t> 77r and the 
bare vertex T 7 tv^a are gauge invariant. However the full amplitude defined by Eq. (^) 
involves off-shell tiN scattering, as defined by Eqs. (|6|)-(|il"|), is not gauge invariant. There 
exists a simple prescription to eliminate this problem phenomenologically. This amounts to 
defining the conserved currents for Eq.(£|) as 



r = J»{SL) - q ' J{SL) n^ (15) 
n ■ q 

where J(SL) is calculated from our model defined above, and n is an arbitrary four vector. 
It is obvious that the currents J M defined by Eq. flip] ) satisfies the gauge invariance condition 
q ■ J = 0. If we use the standard choice of the photon momentum q = (u, 0, 0, | q |) and 
choose n = (0, 0, 0, 1), we then have 

Jo — Jo(SL), 
Jx = Jx(SL), 

Jy = USL), (16) 



and 



J z = J Z {SL) 



ujJo(SL) — 


Q 


USL) 









] ^Jo(SL). (17) 

The above equations mean that within our approach any N(e, eV) observable depending 
on the z-component of the current is determined by Eq. ( ]T7| ) using the time component of 
the SL model, not by J Z (SL). This is very similar to the prescription used in many nuclear 
calculations. We find that the data we have considered in this work can be described by the 
conserved currents J M defined by Eqs. (|T6|) and (|T7|) . We have briefly investigated the model 
dependence due to the freedom in choosing n. We have found that the choice n = (1, 0, 0, 0), 
which leads to J = J{SL) and Jo = ^J Z {SL) gives very similar results at high Q 2 . The 
differences at low Q 2 also appear to be not so large. All results presented in this paper are 
from using the choice Eqs. ([16| ) -([T7| ) . We emphasize that this choice is a phenomenological 
part of our model, simply because we have not implemented any substructure dynamics of 
hadrons into our formulation. In fact this is also the case for all existing approaches based 
on the prescription similar to the form of Eq. ([15]). For example n = q is chosen in the 
recent work by Kamalov and Yang PTJ using the dynamical model developed in Ref. [TH 



There exist other prescriptions to fix the gauge invariance problem, such as those suggested 
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in Refs. [p4|-|36[|. We have not explored those possibilities, since they are also not related 
microscopically to the substructure of the hadrons involved, and are as phenomenological as 
the prescription defined by Eqs.(16)-(17). We will return to this problem when our model 
is further developed, as discussed in section V. 

For our later discussions on the 'jN — > A transition, we define some quantities in terms 
of more commonly used conventions. As discussed in detail in Ref. 0, if we replace the 
ttN propagator G nN in Eqs. ©-CD]) b Y G ^n{ e ) = E-E N {k)-E n (k) witn P denoting the 
principal-value integration, the resulting dressed vertex T^v^a is rea l an d can be directly 
compared with the bare vertex T 7 tv^a- The usual E2/M1 ratio Rem and C2/M1 ratio -Rsm 
for the dressed 'jN — > A vertex are then defined by 



[r^V-^AjMl 



-fiEM — TfTR 1 ' \ L °) 



[f 



p ~ 7 V^AJC72 , 

«sm - 7^ — 1 — • i±yj 

[J- 7 AT^AjAfl 

One can also show that at the resonant energy where the invariant mass W =1236 MeV 
and the tiN phase shift in the P33 channel goes through 90 degrees, the multipole components 
of the dressed vertex F ^v^a are related to the imaginary(imM) parts of the —> tcN 
multipole amplitudes(as defined in Ref. || in the irN P33 channel 

Am = [T^-aWi = Nx lm(M±+), (20) 
Ae = [r« N ^ A ]E2 = Nx Im(£ 3 { 2 ), (21) 
Ac = [r^- A ]c2 = N x Im(5 1 3 { 2 ), (22) 

where N is a kinematic factor. From the above relations, we obtain a very useful relation 
that the E2/M1 ratio Rem and C2/M1 ratio Rsm of the dressed 'jN — ► A transition at 
W = 1236 MeV can be evaluated directly by using the — > nN multipole amplitudes 

_ Im(< 2 ) 

REM ~ i^ 2 ? (23) 

_ Ini(^i{ 2 ) . . 

RSM ~ M^g 5 ) ' (24) 



The formula for calculating the multipole amplitudes and various cross sections from the 
total amplitudes T 77r will not be given here, since they are well documented P-|P2|. 



III. THE RESULTS AT Q 2 = PHOTON POINT 

To determine the 7iV — > A form factors defined by Eq. (|I^D, it is necessary to first fix 
their values at q 2 = by investigating the pion photoproduction reactions. This was done in 
Rcf. [0 by applying the formulation outlined in section II. The first step was to investigate 
the 7rN scattering from threshold to the A excitation region. By fitting the ttN phase shifts, 
the parameters characterizing the strong interaction vertices except the lvNN vertex in Fig. 2 
were determined. The pion photoproduction data were then used to determine Gj\/(0) and 
Ge{0) of the — > A transition(Eq. fll4]) ) and the coupling constant g^NN of the uNN 



vertex of Fig. 2. In Ref. we considered previous pion photoproduction data |42| and the 
LEGS data |IJ of photon asymmetry defined by 

„ da i — daw 



_ (25) 
da± + da\\ ' 

where da±(da\\) are the cross sections with photons linearly polarized in the direction per- 
pendicular (parallel) to the reaction plane. To refine our model, we consider the recent Mainz 
data || here. For the photon asymmetry S, the data from Mainz agree very well with that 
from LEGS. The main improvement we have made is from using the Mainz differential cross 
section(da / dfl) data which are much more precise than the data used in Ref. [[?]. 



With the ttN amplitudes calculated from the Model-L of Ref. [0, we find that the 
Mainz data can be best reproduced by setting Gm(0) = 1.85, Ge(0) = 0.025, and g^NN = 
11.5. These values of Gm(0) and Ge(0) are identical to that determined in Ref. 0. The 
ujNN coupling constant is also only slightly larger than the value 10.5 determined there. 
Such a small change in the determined parameters is due to the fact that the previous 



photoproduction data ||IJ are close to the Mainz data except that their errors are larger. 

Our results for '-/p — > n°p and 7p — > -K + n reactions are compared with the data in 
Figs. 3 and 4 respectively. Clearly the agreement is satisfactory in general. For the 7r° 
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production(Fig. 3), we also show the dependence of the calculated asymmetry X on the 
E2/M1 ratio of the bare jN — ► A vertex. The value REM{bare) = —1.3%, which corresponds 
to Gm(0) = 1.85 and Ge(0) = 0.025, seems to be favored by the data. On the other 
hand, such a dependence is much less for the 7r + production(Fig. 4) since the non-resonant 
interactions play a more important role in this channel. There are still some discrepancies 
with the data. In particular, the differential cross sections at high energies E 1 > 350 MeV 
are underestimated. This could be due to the neglect of the coupling with higher mass 
nucleon resonances and two-pion production channels. For the tt + production(Fig. 4), the 
cross sections at low energies are also underestimated slightly. This could be mainly due to 
the deficiency of our nonresonant amplitude which plays a much more important role in 7r + 
production than in n° production. Possible improvements of our model will be discussed in 
section IV. 

The results presented in the rest of this paper are from calculations with Ga/(0) = 1.85 
and Ge{0) = 0.025. In Fig. 5, we show that the predicted M-f+ 2 and E'i+ amplitudes of 
the 7iV — > nN reactions agree very well with the data from the Mainz98 JZIJ and SM95 
P7| analyses. By using Eqs.([2"3"D-([2"4|) and reading the results at resonant energy E n = 340 
MeV displayed in Fig. 5, we find that the E2/M1 ratio for the dressed ^yN — > A transition 
is Rem = —2.7%. The dotted curves in Fig. 5 are obtained from setting the non-resonant 
interaction v 17T to zero in the calculations. Clearly the non-resonant mechanism has a crucial 
role in determining the electromagnetic excitation of the A. At the resonance energy E 1 = 
340 MeV the bare A amplitude is about 60% of the full amplitude for the Ml transition 
and almost a half for the E2 transition. As discussed in Refs. [0,[np, this large difference 
between the bare A and full amplitudes is the source of the discrepancies between the quark 
model predictions of the 7iV — > A transition and the values determined from the empirical 
amplitude analyses such as those listed by PDG |32|j . 

In Table I, we list our results for the helicity amplitude A3/2 and E2/M1 ratio of the 
syN — > A transition and compare them with the results from other approaches. We see that 
the dressed E2/M1 values from different approaches are very close. For our model and the 
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model of Ref. [pi]] , the large differences between the dressed and bare values are evident, 



indicating the importance of the non-resonant mechanisms in determining the jN — > A 
transition. The bare values are clearly close to the quark model predictions. With Gm{Q) 
and Ge{0) determined, we can then investigate the pion electroproduction reactions. 

IV. PION ELECTROPRODUCTION 

With the matrix element T 77r calculated by using the formula outlined in section II, it is 
straightforward to calculate various observables for pion electroproduction reactions. The 
needed formulation is well documented; see, for example, Refs. [3-|12|. We therefore will 
only give explicit formula which are needed for discussing our results. 

We first consider the unpolarized differential cross sections of the 7*iV — > ttN transition, 



where 7* denotes the virtual photon. In the usual convention pi, it is defined by 



da da T da L r da T da P 

+ + v / 2e(l + e) — - cos fa + e— - cos 2 fa, (26) 



where various differential cross sections depend on the pion scattering angle 8 n , photon 
momentum-square q 2 = —Q 2 = u 2 — q 2 , and the invariant mass W of the final irN system. 
4> n is the off-plane scattering angle between the tt — N plane and the e — e' plane. The 
dependence of Eq.(26) on the angle 9 e between the outgoing and incoming electrons is in 
the parameter e = [1 + ^-i-tan 2 ^Oe]" 1 - Recalling Eq.(2.14) of Ref. 0], we note here that 
the transverse cross section ax and polarization cross section a> are only determined by the 
transverse currents J x and J y and the longitudinal cross section a^ only by the longitudinal 
current J z . On the other hand, the interference cross section 07 is determined by the real part 
of the product J X J*. In general, the contributions from the longitudinal current are much 
weaker than that from the transverse currents. Thus the longitudinal current can be more 
effectively studied by investigating the observables which are sensitive to 07. This has been 
achieved in the recent experiments at MIT-Bates and JLab by utilizing the ^-dependence 
in Eq. (EST). 
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To proceed, we need to define the strength of the charge form factor at q 2 = 0. First 
we use the low momentum limit to set G c (0) = m ^ A G g (0). With G E (0) = +0.025 
determined in section III, we thus have Gc{0) = —0.209. Next we consider the JLab 
p(e, e'rr ) data || at Q 2 = 2.8, 4 (GeV/c) 2 . Since the data are extensive enough, we are 
able to extract each 0^— dependent term in Eq. fl26|). We then adjust Gm(Q 2 ),Ge(Q 2 ) 
and Gc{Q 2 ) to fit the data of these extracted components. Our best fits are the solid 
curves shown in Fig. 6. We also show that the interference cross section d(T//dO(which is 
determined by Re(J x J*)) is sensitive to the charge form factor Gc of the bare — > A 
vertex. In Figs. 7 and 8, we show the comparison with the original data at several off-plane- 
angle <p n . Clearly the agreement is satisfactory. Similar good agreement with the data at 
different W are also found. Some typical results are shown in Fig. 9 for W=1115, 1145, 1175, 
1205 MeV. 

We follow Refs. [0,^1|] to fit the determined form factor values at Q 2 = 0, 2.8 and 4 



(GeV/c) 2 with the following simple parameterization 

G a (Q 2 ) = G a (0)G D (Q 2 )R a (Q 2 ), (27) 

with a = M, E, C and 

Gd{Q2) = (l + Q 2 /0.71(GW/c) 2 ) (28) 
is the usual proton form factor. We find that our results can be fitted by choosing 

R a (Q 2 ) = (1 + aQ 2 ) exp(-6Q 2 ), (29) 

with a = 0.154 and b = 0.166 (GeV/c) 2 for a = M and C. The unpolarized cross section 
data are less sensitive to Ge- Nevertheless, the allowed Ge{Q 2 ) values at Q 2 = 0, 2.8 and 4 
(GeV/c) 2 seem to also follow Eq. (p9"D . For simplicity, we use Eq.(p9"D for a = M, E and C in 
all of the calculations presented below. This simple parameterization is similar to that used 
in Ref. ETJ. Eq. (|27|) then allows us to make predictions for other values of Q 2 . Note that 



with a = 0.154 and b = 0.166 (GeV/c) 2 , the Q 2 -dependence due to R a (Q 2 ) is very small 
compared with the dipole form factor Gd{Q 2 ) in Eq.(p8|). 



13 



To further explore the A excitation within our model, we show in Fig. 10 the Q 2 - 
dependence of the predicted 7*iV — ► ttN multipole amplitudes M^+ 2 , and S^{ 2 at 

W = 1236 MeV where the irN phase shift in P33 channel reaches 90°. Hence, their real 
parts are negligibly small and are omitted in Fig. 10. These amplitudes are proportional to 
the dressed — > A transition strengths A M , A E and A c , as defined by Eqs. (|2Ti|)-([2"2"|). 
We also show the results(dotted curves) from neglecting the non-resonant interaction t> 77r in 
the calculations. It is interesting to note from Fig. 10 that the non-resonant interaction t> 7vr 
enhances strongly these amplitudes at low Q 2 , but much less at high Q 2 . 

From the results(solid curves) shown in Fig. 10 and Eqs. (^) and (p4|), we obtain the 
the Q 2 -dependence of the E2/M1 ratio Rem and C2/M1 ratio R$m for the dressed — > A 
transition. The results are shown in Fig. 11 and Table II. We see that Rsm drops significantly 
with Q 2 and reaches ~ 13% at Q 2 = 4 (GeV/c) 2 , while Rem remains ~ —3% in the entire 
considered Q 2 region. This difference reflects an non-trivial consequence of our dynamical 
treatment of the non- resonant interaction, as seen in Eq.(9). It will be interesting to test 
our predictions in the entire Q 2 region. Clearly, the Q 2 < 4 (GeV/c) 2 region is still far away 
from the perturbative QCD region where Rem is expected to approach unity |3"8fl . 

The dressed 'yN — > A vertex defined by Eq. (|9[) can be cast into the form of Eq. (|TJP at 
the resonance energy W = 1236 MeV. This allows us to extract the dressed Ml form factor 
G* M (Q 2 ) from our results and compare it with the bare form factor Gm(Q 2 ) of Eq. (|14|). 
This is shown in Fig. 12 where we measure their Q 2 -dependence against the proton dipole 
form factor Gd{Q 2 ) defined in Eq.([28|). We see that both bare and dressed Ml form factors 
drop faster than the proton form factor. The difference between the solid and dotted curves 
is due to the non- resonant term(see Eq. (Q)) which can be interpreted as the effect due to 
the pion cloud around the bare quark core. This meson cloud effect accounts for about 40 
% of the dressed form factor at Q 2 = 0, but becomes much weaker at high Q 2 . This implies 
that the future data at higher Q 2 will be more effective in exploring the structure of the 
bare quark core which can be identified with the current hadron models, as discussed in Ref. 
2|. 
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With the 'jN — > A form factors given by Eqs. (p7|)-([29|), we then can test our model 
by comparing our predictions with the data at other values of Q 2 . We first consider the 
MIT-Bates data |§ at Q 2 = 0.126 (GeV/c) 2 . In Fig. 13, we show the results(solid curves) 
for Alt which is defined as(see Eq.(26)) 

= 180°) = 0°) 



tin i ., 1 oAOi tin 



180°) + jg-fa, = 00) 



doT _|_ g do l _|_ gdffp V / 

d^TT d^Tr df27r 

In the same figure, we also show the results(dotted curves) from setting Ge = Gc = for 
the bare 7iV — > A vertex(Eq. (|T^)). The differences between the solid and dotted curves 
indicate the accuracy needed to extract these two quantities of the bare 7iV — > A transition 
within our model. 

Clearly, our predictions are close to the data at W = 1.171, 1.232 GeV. The result at 
W =1.292 GeV appears to disagree with the data. However our model is expected to be 
insufficient at this higher energy, as already seen in the results Figs. 3-4 for photoproduction. 
It is necessary to extend our model to include additional reaction mechanisms such as the 
7iV — > 7rA — > nN transition and higher mass nucleon resonances. 

In Fig. 14, we compare our results(solid curve) with the data for 

q c da T da L 

Rtt = + ^ k = 1800 ' (31) 

where k c and q c are the momenta for the pion and photon in the irN center of mass frame. 
In the same figure, we also show the results(dotted curve) obtained from neglecting the 
nonresonant interaction t> 77r which renormalizes the 'yN — > A vertex and generate non- 
resonant amplitude t 7r7 , as seen in Eqs. (|5|), (||) and (H). The importance of the non-resonant 
interaction is evident. Our predictions reproduce the position of the A peak which is shifted 
significantly by the non-resonant interaction, but underestimate the magnitude at the peak 
by about 15 %. 

In Fig. 15, we present our results for the the induced proton polarization P n for 8^ = 180° 
and the polarization vector n perpendicular to the momentum of the recoiled proton. The 
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data at W =1232 MeV is also from the measurement at MIT-Bates. Our model clearly 
only agrees with the data in sign, but not in magnitude. More experimental data for this 
observable are needed to test the energy-dependence of our predictions. 

The results shown in Figs. 13-15 are from the calculations using the A form factors 
given by Eqs. (0)-© and fitted to the data at Q 2 = 0, 2.8, 4.0 (GeV/c) 2 . The chosen 
parameterization is rather arbitrary and there is no strong reason to believe it should give 
reliable predictions for the A form factors at the considered Q 2 = 0.126 (GeV/c) 2 . We 
therefore have explored whether the discrepancies seen in Figs. 13-15 can be removed by 
adjusting Gm(Q 2 ),Ge(Q 2 ) and Gc(Q 2 ) for the bare 7iV — > A transition. It turns out that 
we are not able to improve our results. It is necessary to also modify the nonresonant 
amplitude t 7vr . We will discuss possible improvements in the next section. 

In Fig. 16, we compare our predictions with some of the Bonn data [fE| at Q 2 = 0.45, 0.75 
GeV/c 2 . Our predictions are in good agreement with these data, but these data have large 
errors. The new high-accuracy data from JLab in this Q 2 region will give a more critical 
test of our predictions. For the forthcoming inclusive p(e, e') data from NIKHEF, we also 
have made the predictions at Q 2 =0.11 (GeV/c) 2 for two polarization observables Att and 
A TL which are clearly defined by Eqs. (2.25b) and (2.25c) of Ref. jn| . Our predictions(solid 
curves) are given in Fig. 17. The dotted curves are from setting the 7*iV — > ttN multipole 
amplitudes E^ 2 and Sf{ 2 to zero. This gives an estimate of the required experimental 
accuracy in using the forthcoming data of A TT and A TL to extract these two amplitudes 
which contain information about Ge and Gq of the 'jN — > A transition. 



V. SUMMARY AND FUTURE DEVELOPMENTS 

In this work, we have extended the dynamical model developed in Ref. |7j to investigate 
th pion electroproduction reactions. The model is first refined at the Q 2 = photon point 
by taking into account the recent pion photoproduction data from Mainz 0. It is found 
that the extracted Ml strength G M (0) = 1.85 and E2 strength G E (0) = 0.025 of the bare 
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7^V — > A vertex are identical to that determined in Ref . |7j . By using the long wave length 
limit, we then obtain Gc{0) = —2.09 for the charge form factor of the bare 7iV — ► A 
transition. 

For the investigation of pion electroproduction, we follow the previous work to define the 
form factor at each photon vertex in the non-resonant interaction t> 77r illustrated in Fig. 2. 
At each value of Q 2 , the — > A transition strengths Gm(Q 2 ), Ge(Q 2 ) and Gc{Q 2 ) are 
the only free parameters in our calculations. We find that the recent p(e, eV°) data at Q 2 = 
2.8, 4 (GeV/c) 2 from JLab @, and at Q 2 = 0.126 (GeV/c) 2 from MIT-Bates [| can be 



described to a very large extent if the bare *yN — > A form factors defined by Eqs. (|27|)-(]29|) 
are used in the calculations. It is found that the remaining discrepancies can not be resolved 
by only adjusting these A form factors. 

We focus on the investigation of the Q 2 -dependence of the A excitation mechanism. It is 
found that the non-resonant interactions can dress the 'jN — > A vertex to enhance strongly 
its strength at low Q 2 , but much less at high Q 2 (Fig. 10). The determined C2/M1 ratio 
(Rsm ) drops significantly with Q 2 and reaches ~ 13% at Q 2 = 4 (GeV/c) 2 , while the E2/M1 
ratio (-Rem ) remains at ~ —3% of the value at the Q 2 = photon point (Fig. 11 and table 
II). The determined Ml form factor drops faster than the usual dipole form factor of the 
proton(Fig. 12). This is in agreement with the previous findings |3"l,|21~ll . 

To end, we turn to discussing possible future developments within our formulation. The 
model we developed in Ref. and applied in this work is defined by the effective Hamiltonian 
Eq. ([3D. It is derived from using a unitary transformation method up to second order in the 
vertex interaction Hj(Eq. (0)). We further assume that the irN and 'yN reactions can be 
described within a subspace A©7r NQ'jN. From the point of view of the general Hamiltonian 
Eqs. ([[])- (0), which can be identified []29| with a hadron model, our model is clearly just a 



starting model. For this reason, no attempt is made here to adjust the parameters of the 
non- resonant interaction t> 77r to perform a x 2 — fit to the electroproduction data. To resolve 
the remaining discrepancies between our predictions and the data, it is necessary to improve 
our model in several directions. 
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To improve our results in the higher energy region(Figs. 3-4,13), we need to include the 
coupling with other reaction channels. An obvious step is to extend the effective Hamiltonian 
Eq. (Q) to include the transitions to i]N, irA and pN states and to include some higher mass 
iV* states. The resulting scattering equations will be more complex than what are given in 
Ref. and outlined in Eqs. fl5|)-(|ITD. In particular, it will have the 7r7riV cut structure due 
to the A — > ttN decay in the ttA channel and the p — > titt decay in the pN channel. This 
7T7tN cut must be treated exactly in any attempt to explore the structure of N* resonances. 
This was well recognized in the early investigations |44| and must be pursued in a dynamical 
approach. 

The second necessary improvement is to use more realistic form factors in defining the 
photon vertices of non-resonant interaction v ln . The prescription Eq. (|T2"D must be relaxed. 
In particular we should use a form factor F ln7T (Q 2 ) predicted from a calculation which 
accounts for the off-mass- shell properties of the exchange pion in Fig. 2. Similarly, the 
vector meson form factor gy ^(Eq. flT3|)) must also be improved since the exchanged vector 
meson is also off its mass-shell. The improvement in this direction could be possible in the 
near future, since the calculations for such off-mass-shell form factors can now be performed 
within some QCD models of light mesons [[EJ . 

This work was supported by U.S. DOE Nuclear Physics Division, Contract No. W-31- 
109-ENG and by Japan Society for the Promotion of Science, Grant-in-Aid for Scientific 
Research (C) 12640273. 
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TABLES 



TABLE I. Helicity amplitude A 3/2 and E2/M1 ratio Rem for the ~/N -> 
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TABLE II. The Q 2 -dependence of the E2/ Ml ratio Rem and C2/M1 ratio R SM for the dressed 



— > A transition calculated from this work. 
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FIGURES 




FIG. 1. Graphical representation of the <-> A interaction. 



/ ' / 




FIG. 2. Graphical representation of the non-resonant interaction of of Eq. (||). The waved, 
dashed, and solid lines denote the photon, pion, and nucleon respectively. 
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£ 7 [MeV] £ r [MeV] 

FIG. 3. The calculated differential cross section(do"/dQ) and photon asymmetry (£) of 
the p(7, 7T°)p reaction are compared with the Mainz data 0. The results from setting 
R EM Q>are) = =0, -1.3, and -2.6 % with G M (0) = 1.85 for the bare -/N -> A vertex are 

indicated in the figure. The three results for the differential cross section are not distinguishable. 
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FIG. 4. Same as Fig. 3 except for the p(j, n + )n reaction. Three results are not distinguishable 
here. 
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FIG. 5. The predicted M l+ and E 1+ amplitudes for the jN — > ttN reaction are compared 
with the results from the empirical amplitude analyses. The dotted curves are from the calculations 
neglecting the non-resonant interaction text for more detailed description. The open circle 



data are from SM95 37] and solid circle data are from Mainz98 [26| 
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FIG. 6. Three components of the calculated p(e,e'7r°) differential cross section(Eq. at 
Q 2 = 2.8(upper row), 4.0(lower row ) (GeV/c) 2 and W = 1235 MeV are compared with the data 
which are extracted from the JLab data ||(some of these original data are shown in Figs. 7-9 ). 
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FIG. 7. The predicted p(e,e'7r°) differential cross sections at W = 1235 MeV, Q' 2 
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(GeV/c) 2 are compared with the JLab data 
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FIG. 9. The predicted p(e, e'ir°) differential cross sections at <p n = 135° with Q 2 = 2.8(left) and 



4(right) (GeV/c) 2 are compared with the JLab data |J. 
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FIG. 10. Q 2 dependence of the imaginary(Im) parts of the 7*iV — ► irN multipole amplitudes 
M x 3 { 2 E^l 2 , and at W = 1236 MeV. The solid curves are from our full calculations, and the 
dotted curves are from calculations with the non-resonant interaction set to zero. The real 
parts at W=1236 MeV are negligibly small and are omitted. 
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FIG. 11. The predicted Q 2 -dependence of the E2/M1 ratio i? EM and C2/M1 ratio R SM of the 
dressed 7iV — ► A form factors. 
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FIG. 12. The ratio between the Ml form factor of the jN — ► A transition and the proton 
dipole form factor Gd defined by Eq. (|28|). The solid curve is G* m (Q 2 )/Gd(Q 2 ) for the dressed 
Ml form factor, and the dotted curve is Gm(Q 2 )/Gd(Q 2 ) for the bare Ml form factor 
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FIG. 13. The predicted asymmetry Alt (Eq. (p0|) ) of the p(e,e'7r°)p reaction at 
Q 2 = 0.126(GeF/c) 2 are compared with the data from MIT-Bates ||. The dotted curves are 
obtained from setting Gm{Q 2 ) = Gc{Q 2 ) = 0. 
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FIG. 14. The predicted Rtt(E>Q- (31) at 6 W = 180° are compared with the data from 
MIT-Bates pi. The dotted curve is obtained from setting the non-resonant interaction Vj n to zero. 
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FIG. 15. The predicted induced proton polarization P n at 9 n = 180° and the polarization vector 
n perpendicular to the recoiled proton momentum are compared with the data from MIT-Bates 
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FIG. 16. The predicted differential cross sections of p(e,e'n°) reaction at Q 2 =0.45(left), 
0.75(right) (GeV/c) 2 and W =1232 MeV are compared with the data @. 
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FIG. 17. The predicted A TT and A TL , as denned by Eqs. (2.25b) and (2.25c) of Ref. 0, for 
the inclusive p(e,e') reaction at Q 2 = 0.11 (GeV/c) 2 . The dotted curves are obtained when the 
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7 irN multipole amplitudes S l+ and EJ + s are not included in the calculation. 
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